Magneto-elastic coupling between copper spin configurations and oxygen octahedra in 
La2-a;Sr2:Cu04. Untwinning, Raman (phonon) spectrum, and neutron response 
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We reinterpret neutron scattering and Raman spectroscopy experiments in La2-a;Sr2!Cu04, show- 
ing modulation of spin and lattice degrees of freedom both in the low-temperature orthorhombic and 
tetragonal phases, in terms of a magneto-elastic coupling between noncollinear configurations for 
the copper spins and the oxygen octahedra. The magneto-elastic coupling furthermore explains the 
recently discovered magnetic-field induced untwinning in La2-a:Sr2,Cu04, for x < 2%, and allows 
us to understand why in ha.2-x-y^dySTxC'a04 with y = 0.4 and x = 1/8, for example, when (1,0) 
and/or (0, 1) static spin spirals are stabilized due to hole frustration, the atomic modulation in the 
low-temperature tetragonal phase of La2-x-t/Ndi,Sra;Cu04 has a periodicity of 4 lattice spacings, 
half the magnetic one. 
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PACS numbers: 74.25.Ha, 75.10.Jm, 74.72.Dn 



INTRODUCTION 



La2-xSra;Cu04 is one of the best studied high tem- 
perature superconductors and is a system where the 
phenomenon of ID charge segregation, or stripes for- 
mation, has been argued to occupi. In fact, inelas- 
tic neutron scattering experiments within the supercon- 
ducting phase of La2-a;Srj;Cu04, x > 5.5%, have re- 
vealed that dynamical incommensurate (IC) spin corre- 
lations coexist with superconductivitji^. In addition, it 
has also been observed experimentally the existence of 
IC "charge" peaks accompanying the IC magnetic or- 
der (even though neutrons cannot directly probe charge 
distributions) , with twice the incommensurability of the 
magnetic one. This phenomenon is most clearly observed 
when a low-temperature tetragonal (LTT) phase is sta- 
bilized, for example by Nd doping^ and it was imme- 
diately interpreted as an evidence of the formation of 
charge stripes that act as antiphase domain walls. ^ In the 
same spirit, the observation of diagonal static IC mag- 
netic correlations at even lower doping, 2% < a; < 5.5%, 
within the low-temperature orthorhombic (LTO) spin- 
glass (SG) phase of La2-a;Sra;Cu04f4 as well as recent 
evidence of diagonal "charge" modulation from Raman 
spectroscopyjSi have also been used to support the pic- 
ture of diagonal stripe formation at lower doping. 

In this article we propose an alternative interpretation 
for the above described experimental results as originat- 
ing from a magneto-elastic coupling between the Cu"'""'" 
spins and the lattice in La2-2:Sr2:Cu04.^ We start by ex- 
plaining the recently discovered magnetic-field induced 
untwinning in La2-2;Sra;Cu04ii According to Ref. 0, 
the two inequivalent orthorhombic directions of La2Cu04 
can be swaped at room temperature by the application 
of a strong in-plane magnetic field, _B ~ 14 T. This re- 
markable effect was confirmed by X-ray measurements 
showing changes in the structure of La2-a;Srj;Cu04, and 
the crystallographic rearrangement observed in an optical 
microscopeii This effect has been used since then to pro- 



duce extremely high quality untwinned La2-a;Srj;Cu04 
single crystals. We then argue that it is precisely this 
same magneto-elastic coupling that gives rise to a super- 
modulation of the oxygen octahedra when a noncollinear 
configuration for the Cu^"'' magnetic moments is stabi- 
lized, as a consequence of the magnetic frustration in- 
troduced by hole-doping)^ and we discuss its effects on 
the Raman and neutron responses, for the case of non- 
collinear magnetic structures both within the LTO and 
LTT phases of La2-xSr^Cu04. 



II. THE MAGNETO-ELASTIC COUPLING 

In La2Cu04 each Cu++ ion is located at the center 
of an oxygen octahedron, see Fig. The crystal field 
produced by the ions of the octahedron splits the 
energy levels of the Cu"'"'" orbitals in such a way that the 
unpaired electron is found at the dx2_y2 orbital. This 
is a planar orbital located at the base of the octahedron 
(the common base of its two pyramids). The spin-orbit 
(SO) interaction in La2Cu04 then introduces a large XY 
anisotropy that favors a configuration for the Cu"'""'" spins 
parallel to the plane of the d^2_y2 orbital. In case there is 
no tilting of the octahedra, this corresponds to confining 
the spins to the Cu02 layers, as it happens for exam- 
ple in the Sr2Cu02Cl2 system. On the other hand, for 
La2Cu04 the tilting of the oxygen octahedra in the LTO 
phase generates a Dzyaloshinskii-Moriya interaction that 
determines an easy-axis for the antiferromagnetic order- 
ing and leads to the canting of the spins out of the Cu02 
planes, see Fig.^ This can be better understood by us- 
ing the notation introduced in Ref. 0. We first write the 
spin of Cu^"'" in terms of its staggered, n, and uniform, L, 
components, S^/S* — e''^ '''n(xi)-f L(xi), where, as usual, 
Q = (tt, tt) is the antiferromagnetic ordering wave vec- 
tor. Then we introduce a thermodynamic Dzyaloshinskii- 
Moriya vector D+ pointing along the axis of tilting of the 
oxygen octahedra, see Fig. 
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The magnetic free-energy density for uniform config- 
urations of the staggered moments in a single layer of 
La2Cu04 is^ 

J^Af = TT- {(n ■ + TxY ni + {n- Bf - 2B ■ (n x D+)} , 

(1) 

where gc — SJa^, Txy is the XY anisotropy energy, J 
is the in-plane antiferromagnetic superexchange, and a 
is the lattice spacing. We now introduce a lattice free- 
energy density 

J^L = -C D+ • iC, (2) 

which favours D-|_ oriented along a certain direction X 
in a fixed laboratory reference frame. The parameter 
C introduced here controlls the softness of the oxygen 
octahedra and in the case C ^ oo the octahedra becomes 
completely rigid and fixed to the laboratory frame. 

Since the orientation of D+ defines the local or- 
thorhombic coordinate system of the crystal, we conclude 
that, at zero magnetic field, B = 0, the equilibrium con- 
figuration of the system is such that D-|_ || a || X with 



n(x,) • D+ = 0, and L(x,) = — (D+ x n(x,)) . (3) 

Together with the large XY anisotropy, which confines 
n(xi) to the orthorhombic ab basal-plane, the above two 
conditions uniquely determine the easy-axis and canting 
of the Cu"'""'" spins, with n || b and L || c, see Fig.Q] 



LTO LTT 




FIG. 1: (Color online) Left: LTO (abc) and LTT {xyz) coor- 
dinate systems. Center and Right: Tilting of the octahedron 
in the LTO (center) and LTT (right) phases. Uniform spin 
component L (blue arrow), staggered spin component n (red 
arrow), and Dzyaloshinskii-Moriya vector D+ (green arrow). 



III. MAGNETIC UNTWINNING 

The physics of the magneto-elastic coupling in 
La2-3;Sr3;Cu04 is essentially given by Eqs. Q and 
These equations imply that the vectors n, L, and D_|., 
form a triad of vectors mutually orthogonal that can only 
move rigidly, see Fig. Let us now consider twins A 
and B such that the orthorhombic basis-vectors in twin 
A are rotated by 7r/2 relative to the ones in twin B, 



'X-A -L Xb, see Fig. El When a strong in-plane mag- 
netic field, B > 2gcQ, is applied parallel to the X^ or- 
thorhombic direction of twin A, we find that the total free 
energy Tm + is minimized, within such twin domain, 
for a new configuration of n and D-|. such that n || ca 
and D_|- || bA (see Fig. [2l middle pannel). This is only 
possible if accompanied by a 7r/2 rotation of the entire 
octahedra in twin A (remember that D_|_ determines the 
axis of the tilting of the octahedra), which then causes 
the swap of the local orthorhombic directions. In twin B 
the effect of the same magnetic field is harmless, since it 
is longitudinal, B _L X^. In fact, a longitudinal magnetic 
field in La2Cu04 causes a continuous rotation of the two 
vectors, n and L, within the be plane (ifliii as a result 
of the term — 2B • (n x D+) in Eq. Since in twin 
B the D+ vector remains perpendicular to the be plane, 
such rotation can take place without causing structural 
changes. When the field is removed, n and L return to 
their original positions (as originally in twin B), but the 
orthorhombic axis in twin A remain swaped. 



B = B = 14T B = 




FIG. 2: (Color online) i) B = (left): the two twins have 
basis-vectors rotated by 7r/2 relatively; ii) B = 14 T: the 
basis-vectors in twin A swap due to the flop of staggered fleld 
n and orient identically to the ones of twin B, where the 
staggered n and uniform L vectors are rotated in the be plane 
(see discussion in the text); iii) B = (right): the staggered 
n and uniform L vectors return to their equilibrium positions 
and the crystal remains untwined. 



The important message from this magneto-elastic al- 
lowed untwining in La2-a;Sra;Cu04 is that a change in the 
spin configuration can have a strong effect on the equi- 
librium position of the tilting axis of the oxygen octahe- 
dra. What is most surprising, however, are the associated 
energy scales, especially being an effect originated from 
SO coupling. Besides the ability of causing the swap of 
the crystallographic axis in a magnetic field (the sam- 
ple actually expanded about 1% along the orthorhombic 
b direction after the untwinningZ') magnetic susceptibil- 
ity experiments in La2-a;Sra;Cu04 by Lavrov et al.^ at 
much lower fields, have shown that the Cu"'""'" spins re- 
main confined to the orthorhombic (6c) plane, for differ- 
ent doping concentrations, < a; < 4%, and for a rather 
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wide temperature range, < T < 400 Km^ This is an- 
other consequence of Eqs. ^ that dictates that the n and 
D+ vectors should be perpendicular A In what follows we 
shall discuss the effects of the magneto-elastic coupling 
on the oxygen octahedra and its effect on the Raman and 
neutron responses, for the case of noncollinear magnetic 
structures both within the LTO and LTT phases. 

IV. RAMAN SPECTRUM IN THE LTO PHASE 

In a recent publication^i^ it was shown that the mag- 
netic frustration introduced by hole-doping in the LTO 
phase of La2-a;Sra;Cu04 gives rise, above a certain criti- 
cal doping, to a static helicoidal supermodulation where 
the n vector precesses (in space) around the b easy-axis 

n(xi) = (cTa cos (Qs • Xi), nb{Ki), dc sin (Q5 • x^)), (4) 

with the wave vector for the spin incommensurability 
Qs II ^, < CTa < 1, and ri^ = \-n\ - n^y^ Such heli- 
coidal structure was then shown to be consistent with the 
magnetic susceptibility experiments in La2-a;Sra;Cu04 
for X = 2%, 3% and 4%, since it guarantees that the 
total magnetization in a applied field remains confined 
to the he planei^ 

The physical picture behind the formation of such he- 
licoidal structure can be easily understood as follows. 
Above a certain critical doping, that is essentially deter- 
mined by the strength of the Dzyaloshinskii-Moriya gap, 
/^dm/J ~ 0.02, the magnetic frustration introduced by 
the doped-holes favours a noncollinear configuration for 
the Cu"'""'" spins. Since these are coupled to the oxygen 
octahedra by the magneto-elastic coupHng in Eq. lO, 
any spin rotation should be accompanied by a rotation 
of the tilting axis of the octahedra (D+ is now D_|_(xi)). 
However, in the LTO phase of La2-2;Sr2;Cu04, there are 
no n-fold axis of rotation with n > 2, and thus the or- 
thorhombicity frustrates any such rotation. For example, 
the in-plane spiraling of the n vector is forbidden because 
it would require a continuous rotation of the tilting axis 
of the octahedra around the c axis and this is not al- 
lowed in the LTO phase. Thus, for Sr concentration well 
above the critical doping, deep inside the SG phase, the 
resulting noncoUinear spin structure in the LTO phase 
corresponds to, for example, a helicoidal modulation of 
the n vector around the b easy-axis, as in Eq. 1^, ac- 
companied by a precession of the octahedron around its 
equilibrium (x — 0) position, see Fig. El Furthermore, 
such atomic modulation has wave vector for atomic in- 
commensurability given by Qc || b. 

It is worth pointing out that, for even lower doping, 
X « 0.01, 0.024, it has been recently shown by Luscher et 
al^ that the appropriate noncoUinear spin structure is 
such that the n field instead oscillates within the CUO2 
plane around the b easy-axis. We emphasize here that, 
once the magneto-elastic coupHng in Eq. (pj is taken into 
account, such oscillation also gives rise to a modulation 
of the octahedra, just as in the case of the helicoid. 




FIG. 3: (Color online) Precession of the triad (left) and oc- 
tahedron (center) around their equilibrium {x = 0) position. 
The precession takes place as one goes along the b axis, such 
that the two wave vectors Qs and Qc are parallel to b. Right: 
neutron response in the LTO phase due to modulation of n. 
The neutron response is shifted from the N6el peak (unfilled 
square at (7r,7r)) to new diagonal incommensurate positions 
(filled circles), due to the helicoidal magnetic modulation de- 
scribed by Eq. 0] 



The magnetic modulation with wave vector Qs || b, 
either helicoidaU^ or oscillatorjiii, gives rise to the diag- 
onal IC peaks observed in the neutron response^ see Fig. 
El A supermodulation of the octahedron, in turn, would 
cause the enlargement of the atomic unit cell. This would 
result in the appearance of a plethora of extra vibrational 
modes that could be accessed via Raman spectroscopypi^ 
As we shall now demonstrate, this is actually the case. 
On the top of Fig. 31 x = 0, we observe 5 sharp phonon 
peaks, corresponding to the symmetric vibrational modes 
predicted by group theory and accessible in the (AA) 
(red-solid-curve) and (BB) (blue-dashed-curve) Raman 
geometries. We also observe a sharp Big phonon signal 
at 220 cm~^, accessible in the {BA) (black-dotted-curve) 
geometry. As doping concentration is increased, x = 1% 
and 3%, the 5 sharp phonon peaks are replaced by much 
more numerous and broader ones while the Big phonon 
remains sharp, although with decreasing intensity, see 
Fig. 01 This is a consequence of the enlargement of the 
atomic unit cell associated with the supermodulation of 
the octahedra with wave vector Qc || b (Qc is the wave 
vector for the modulation of D+(xi)), see Fig.O 

In fact, recent Raman experiments in La2-a;Sra;Cu04 
at X = 2% showed indications of anisotropic atomic mod- 
ulation in the B2g Raman geometryi^ These results are 
consistent with the above physical picture of the octa- 
hedra supermodulation. The atomic modulation result- 
ing from the precession of the octahedra has wave vector 
Qc II b. In order to probe such modulation, we must 
choose a Raman geometry such that the wave vector Qc 
has a finite projection on the axis of the electric field of 
both incoming em and outgoing Cout Hght. Thus, the 
above octahedra supermodulation should become acces- 
sible exactly in the i?2g geometry, where Cm and Cout are 
oriented along the Cu-Cu bonds, with respect to which 
Qc is diagonal.^ Furthermore, it should not be acces- 
sible at all in the (BA) geometry (black-dotted-curve) 
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FIG. 4: (Color online) Raman response at 10 K in the {AA) 
(red-solid), (BB) (blue-dashed), and (BA) (black-dotted) po- 
larization configurations, for a; = 0, 1%, and 3% (from Ref. 
IT^ . For X = we clearly see sharp phonon peaks corre- 
sponding to the 5 symmetric Raman active modes predicted 
by group theory, as well as one Big phonon. For x = 1% and 
3%, however, a larger number of peaks is observed, signaling 
the enlargement of the atomic unit cell. 



where Cm and Cout are oriented along the orthorhombic 
basis- vectors, and this results in a rather featureless Big 
response (except for the phonon at 220 cm~^), see Fig.^ 
Although the results by Tassini et al.^ were interpreted as 
a signature of diagonal stripe formation, we beheve that 
the octahedra supermodulation provides a more reason- 
able explanation for the Raman response reported in Fig. 
01 and in Ref. 0, especially for the x = 1%,3% samples 
where stripes can be safely ruled out. 



V. NEUTRON RESPONSE IN THE LTT PHASE 

When a tetragonal phase is stabilized in 
La2-a;Sra;Cu04, by Nd doping for example, incom- 
mensurate "charge" peaks are observed in neutron 
scattering, whose incommensurability is twice the 
magnetic one, see Fig. Even though neutrons cannot 
directly couple to charge excitations, this result has been 



interpreted, for over 10 years, as a sign of charge stripe 
formation acting as anti-phase domain walls. 

Let us now demonstrate how the magneto-elastic Eqs. 
(PJ can be used to explain the neutron scattering results 
within the LTT phase. In the LTT phase, the D+ vector 
points along one of the basal tetragonal directions, say 
X, such that two oxygen atoms remain confined to the 
tetragonal basal plane while the other two are tilted, one 
above and one below the Cu02 layers, see Fig. ^ From 
the magneto-elastic coupling equations Eq. (EJ we then 
find n II y and L || z. There is one major difference, 
however, regarding possibilities for the rotation of the 
tilting axis of the octahedron between the LTO and LTT 
phases when a noncollinear spin structure is favoured. In 
the LTT phase, the z axis is in fact a n-fold axis with 
n > 2, and a in-plane rotation of the staggered moments 
is now not only allowed but energetically favourable (as 
discussed earlier, the XY anisotropy favours n(xi) to be 
parallel to the planes) . The above argument is consistent 
with the calculations by Sushkov and Kotovjl- using real- 
istic parameters for the i— t'— t" — J model, which showed 
that the magnetic ground state, in the LTT phase, is 
composed by (1,0) and/or (0,1) spirals 



n(xi) = (cos (Qs • Xi), sin (Qs • Xi), 0), 



(5) 



with Qs ^ 6 ^ X. The magneto-elastic coupled spiraling 
of the Cu''"'" spins and tilting axis of the octahedra, D+, 
produces the pattern depicted in Fig. [3 

Now, by analyzing Fig. [^we can finally understand 
why the atomic incommensurability is twice the magnetic 
one, and why for x = 1/8 the atomic periodicity is of 4 
lattice spacings. Since Qs = S ^ x, for x — 1/8 the 
spin periodicity is of 8 lattice spacings, obtained after 
8 rotations of 7r/4pi^ However, after 4 rotations of tt/A, 
when the spins are anti-parallel, the configuration of the 
octahedron is equivalent to the initial one, see Fig. 
Although this might not be apparent on a first inspection 
due to the difference in the tilting of the oxygen atoms, 
we should not forget that in a body-centered tetragonal 
unit cell, such as the one of La2-a;Sra;Cu04, corner and 
central Cu"*'"'' ions are equivalent. Thus, any symmetry 
operation can be defined together with a half translation 
of the lattice vectors. This brings the octahedron on the 
far right of Fig. to the one in the far left, and so they 
are indeed equivalent. Notice now that the modulation of 
the atoms in the octahedra is half the magnetic one, and 
so the neutron response in the reciprocal space should 
be such that the atomic incommensurability is twice the 
magnetic onep see Fig. 

Finally, the magneto-elastic coupling might also gives 
us a hint towards the understanding of some un- 
expected features on the phonon spectrum of the 



La2 



^NdySra;Cu04 system, at x = 1/8 and y = 0.4, 



recently reported by the Raman experiments of Gozar 
et al. in Ref. It has been observed that the lowest 
energy of the Ag phonon modes, there denoted A mode, 
has a rather broad Raman spectral line. Since this is 
the phonon mode related to the spatial distribution of 
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Atomic supermodulation in the LTT phase 
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FIG. 5: (Color online) Top: Magneto-elastic coupled spiral- 
ing of the Cu"*"^ spins (black arrows, V-shaped tip) and D+ 
(green arrow, triangular tip). Blue (larger and hatched) oxy- 
gen atoms are canted up the plane, red (smaller and filled) 
ones are canted down the plane, and open ones are confined 
to the basal plane. The octahedra on the far left and far right, 
exactly at half the magnetic periodicity, are equivalent (see 
discussion in the text). Bottom: Since the atomic periodicity 
is half the magnetic one, the neutron response in reciprocal 
space is such that the atomic incommensurability (large filled 
squares) is twice the magnetic one (smaller filled circles). 



the tilting angle of the oxygen octahedra, the Raman 
data seems to indicate a locally fluctuating octahedra tilt 
distribution in La2-x-y^dySTxCu04r^ that can be un- 
derstood in terms of short-range magneto-elastic coupled 
spin-tilt-axis fluctuations. 



VI. CONCLUSIONS 

We have shown that neutron scattering and Raman 
spectroscopy experiments, either in the LTO or LTT 
phases of La2-3;Sra;Cu04, as well as the process of un- 
twining in a strong in-plane magnetic fleld, can be un- 
derstood in terms of the unusually strong magneto-elastic 
coupHng between the noncoUinear conflguration for the 
Cu"^+ spin and the oxygen octahedra. Although the 
analysis here perfomed can go no further than discussing 
static phenomena of charge and spin ordering, it is very 
much possible that the distribution of oxygen doped holes 
in La2-2;Sr2;Cu04 will follow the anisotropic ID octahe- 
dra modulation proposed above, thus presenting us with 
an interesting possibility for the origin of the anisotropic 
ID electronic (charge) dynamics in cuprates. This inter- 
esting possibility is presently under investigationiii 
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